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ABSTRACT 

The newly developed X-ray visibility forward fitting technique is applied to Reuven Ramaty High 
Energy Solar Spectroscopic Imager (RHESSI) data of a limb flare to investigate the energy and height 
dependence on sizes, shapes, and position of hard X-ray chromosphcric footpoint sources. This pro- 
vides information about the electron transport and chromospheric density structure. The spatial 
distribution of two footpoint X-ray sources is analyzed using PIXON, Maximum Entropy Method, 
CLEAN and visibility forward fit algorithms at nonthermal energies from ~ 20 to ~ 200 keV. We 
report, for the first time, the vertical extents and widths of hard X-ray chromospheric sources mea- 
sured as a function of energy for a limb event. Our observations suggest that both the vertical and 
horizontal sizes of footpoints are decreasing with energy. Higher energy emission originates progres- 
sively deeper in the chromosphere consistent with downward flare accelerated streaming electrons. 
The ellipticity of the footpoints grows with energy from ~ 0.5 at ~ 20 keV to ~ 0.9 at ~ 150 keV. 
The positions of X-ray emission are in agreement with an exponential density profile of scale height 
~ 150 km. The characteristic size of the hard X-ray footpoint source along the limb is decreasing with 
energy suggesting a converging magnetic field in the footpoint. The vertical sizes of X-ray sources 
arc inconsistent with simple collisional transport in a single density scale height but can be explained 
using a multi-threaded density structure in the chromosphere. 
Subject headings: Sun: flares - Sun: X-rays, gamma rays - Sun: activity -Sun: particle emission 



1. INTRODUCTION 

In the standard flare scenario, electrons that were 
accelerated in the corona stream downwards toward 
the dense layers of the solar atmosphere, where they 
are stopped via collisions producing intense hard X- 
ray (HXR) emission in the chromosphere. Higher en- 
ergy electrons penetrate deeper into the chromosphere. 
Therefore, measurements of the spatial structure of the 
HXR emission as a function of energy provide informa- 
tion about the chromospheric density. Being optically 
thin, X-rays give the most direct information about the 
spatial and energy distribution of energetic electrons in 
the solar atmosphere. Prior to the launch of the Ra- 
maty High Ener gy Solar Spectroscopic Imager (RHESSI) 
(jLin et al.ll2002D . HXR instruments typically had l imite d 
imaging-spectroscopy capabilities (jKosugi et al.l I1992D . 
RHESSI's aptitude to image in various energy ranges 
opens new horizons for studying the detailed structure 
of HXR emitting sources in the chromosphere. 

RHESSI docs not directly image the Sun but uses 
9 pairs of Rotating Modulation Collimators (RMCs) 
to time-modulate spatial informat ion in the signal ob - 
tained in its germanium detectors (|Hurford et al.|[2~002[ ). 
Each RMC has a different thickness of its slits and 
slats making it sensitive to different spatial scales, pro- 
viding modulation at nine spatial frequencies. The 
reconstruction of an image from these time modu- 
lated lightcurves, can be accomplished by various imag- 
ing a l gorithms (lEmslie et al.l 12 003: B attaglia fc Benj 
20071 IKrucker fc Lid 120081 ISaint-Hilaire et al.1 12001 
Dennis fc P ernak 2009). The new visibility based ap- 



and the phases in the individual roll bins are X-ray visi- 
bilities. This effectively provides two dimensional spatial 
Fourier co mponents (X-ray visibilities) over a wide r ange 
of energies (iHurford et al.(2002HSchmahl et al.l200l . To 
convert the time-modulated signal or X-ray visibilities to 
an image is an in verse problem (e.g. IPiana et al.l 120071 : 
iPrato et al.1 120091 ). The reconstructed images face un- 
avoidable difficulties due to measurement errors, finite 
coverage in Fourier space, and ill-posedness of the recon- 
struction problem. The resulting reconstruction errors 
and small dynamic range makes it difficult to accurately 
measure source sizes from reconstructed images. At best 
the imaging resolution is down to 2 arcseconds but in 
practice is around 7 arcseconds for the typical flare non- 
thermal energy range. 

However, the moments of X-ray source distribu- 
tion (source position, source size, etc) can be inferred 
with hi gher precision either from the time-modulated 
signal (lAschwanden et al.l [200l IKrucker fe Lid" [200l 
ISaint-Hilaire et al.ll2008tlFivian et al.)l2009D o r visibilities 



proach to RHESSI imaging starts by summing (stacking) 
the lightcur ves per roll bins over a few spin periods of the 
spacecraft (jSchmahl et al.l 120071 ). The fitted amplitudes 



(IXu et all I2008t i Kontar et al.l 12008: D ennis fc Pernakl 
I2009t IPrato et al.ll2009t) . Thus. RHESSI measurements 
of X-ray source positions can recover sub-arcsecond 
information using RH ESSI modulated lightcurves 
(jAschwanden et all 120021: ILiu et all [20 06: Mrozek 2006) 
or visibilitie s (IKontar et al.ll2 008: De nnis fc Pernakll2009t 
IPrato et al.l l2009f ). This has allowed clear demon- 
stration of the height-energy dependence of HXR 
sources: higher energy s o urces originate at lower 
heights (lAschwanden et al.l 120021: IBrown et ahl 120021: 
ILiu et all 120061 ). This has substantiall y improved upon 
previ ous results with Yohkoh/HXT (Matsushit a et al.l 
I1992D. The recently d eveloped v isibility-based techni que 
(jSchmahl efail I2007D . allowed IKontar et ahl (J2008D to 
improve previous measurements and infer characteris- 
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tic sizes (FWHM) of HXR footpoints and not just the 
centroid height. From this the convergence of the mag- 
netic flux and neutral hydrogen density distribution in 
the chromosphere was inferred. 

In this paper we study the structure of solar flare hard 
X-ray sources using various imaging algorithms: PIXON, 
MEM, CLEAN and visibility forward fit. Using RHESSI 
X-ray visibilities we find the characteristic shapes and 
positions for different energy ranges. We show that the 
technique of forward fitting X-ray visibilities allows us to 
determine not only the FWHM of the sources but vertical 
and horizontal sizes of the sources, which is required for 
examining the density structures of the chromosphere. 
Theoretical relationships were compared with observa- 
tions to find the density structure of the chromosphere. 
The vertical size of the X-ray sources is found to be larger 
than the ones predicted by a hydrostatic atmosphere in 
thick-target scenario. However, assuming that the elec- 
trons are propagating along several narrow threads with 
different density profiles can explain the measured verti- 
cal sizes of the sources. 
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Fig. 1. — Lightcurves of the 6- Jan-2004 flare. The vertical dotted 
lines show the accumulation time interval 06:22:20-06:23:00 UT 
which is used in the subsequent spectral and imaging analysis. 



rors to forward fit parameters of the model in Equation 
((2| . Reliable error estimates for images reconstructed 
with other algorithms are currently unavailable. 



2. X-RAY VISIBILITIES AND CHARACTERISTIC 

SIZES 

The spatial information about an X-ray source mea- 
sured by RHESSI for a given energy range and time inter- 
val ca n be presented (|Hurford et al.l l2002; Sch mahl et al.l 
12007ft as two dimensional Fourier components or X-ray 
visibilities 



V(u,v;e) 



I(x,y;e)e 27ei ^ u+yv) dxdy (1) 



where I(x, y;e) is the observed image at photon energy 
e. Then, reconstructed X-ray image I(x, y; e) is the in- 
verse Fourier transformation of measured X-ray visibili- 
ties V(u, v; e). Each of the nine RHESSI Rotating Mod- 
ulating Collimators (RMC) measures V(u, v; e) at a fixed 
spatial frequency (or a circle in the (it, v) plane) corre- 
sponding to its angular resolution and with a position 
angle along the circles, which varies continuously as the 
spacecraft rotates. Nine detector grids with angular reso- 
lutions growing with detector number are logarithmically 
spaced in the (u,v) plane. Since the measured visibili- 
ties sparsely populate the (u,v) plane and have statis- 
tical uncertain ties, the direct inve r se Fourier transform 
is impractical jHurford et al.ll2002t iSchmahl et~aTll2007t 
iMassone et al.ll2009ft and alternative methods should be 



used. 

Assuming a characteristic shape of X-ray source, one 
can find the position and characteristic sizes directly by 
fitting a 2D Fourier image of the model to the RHESSI 
visibilities. Here, we assume that the sources can be 
presented as elliptical Gaussian sources 
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where 2V21n2cr x and 2^2 In 2a y are FWHMs of an el- 
liptical Gaussian source in x and y direction respectively, 
Xo(e), 2/o(e) is the position of the source, and Iq is the to- 
tal photon flux of the source. One major advantage of the 
visibility forward fit approach is that knowing the errors 
on visibilities V(u, v; e) one can readily propagate the er- 



2.1. The shape of footpoints for a limb event 

Using HXR data from RHESSI we analysed a limb 
event on January 6t h, 2004 (G OES M5.8 class). As 
shown previously by (jKontar et al. 2008]), this event is 
ideally suited for our analysis having two well sepa- 
rated footpoints: one bright and a second much weaker 
footpoint. In addition, the locati on of the flare at the 
limb greatly reduces a lbedo flux (|Bai fc Ramatvl 119781 ; 
iKontar fc Jeffrev| | 2010f) . so that the albedo correction 
( Kontar e t al. 2006]) becomes negligible. 

The flare occurred at the eastern limb near (—975, 75) 
arcseconds from the disk center at ~ 06 : 22 UT (Fig- 
ure [1]). It was imaged during the time of peak emission 
> 50 keV 06:22:20-06:23:00 UT indicated by the vertical 
dotted lines in Figure [TJ usin g four different ima ge algo- 
rithms (see Figure^): Clean (IH urford et al.ll2002T). MEM 
NJIT (ISchmahl et all 12007ft . PIXON (IPin a fc Puetterl 



119931: iMetcalfet all 11996ft and vi sibility forward fit 
(jHurford et all 12002; ISchmahl et alj|200l . The result- 



ing images in five energy bands covering the nonthermal 
emission are shown in Figure [21 Each image was made 
using the front segments of detectors 2 to 7. Grid 1 
with the highest spatial resolution had no significant sig- 
nal and grids 8-9 are too coarse for our flaring region. 
Previously the flare was imaged using ten energy bins 
(jKontar et al.ll2008ft and simple circular gaussian fit but 
this was reduced to five wider bins in this paper to im- 
prove signal to noise. Figures 12131 demonstrate that the 
brighter source has an elliptical shape at various ener- 
gies, so an elliptical Gaussian could be used as natural 
X-ray distribution model (Figure [3]). 

Comparing the different algorithm results we find that 
CLEANed images have systematically larger sizes than 
the other algorithms. This is related to the fact that 
CLEAN images are determined by the user choices for 
analysis (clean beam size) and not the requirements of 
the data and hence should be used with great care to 
measure source sizes. Q MEM NJIT has produced 

1 Reduction of the CLEAN beam size by 1.7 produces images 
with spatial characteristics similar to other algorithms. The cur- 
rent version of clean does not have a robust procedure to determine 



The sub-arcsccond hard X-ray structure of loop footpoints in a solar flare 



90 p 

s 



• 









- 1 000 -980 96C 1 JOll 98Q 930 



1 000 -98C 



Fig. 2.— HXR images of the January 6th 2004 flare (accumulation time interval 06:22:20-06:23:00 UT) in 5 energy ranges: 18-25, 25-29, 
29-43, 43-75, 75-250 keV (energy grows from left to right) using different algorithms: (from top row to bottom) CLEAN, PIXON, MEM- 
NJIT, VIS forward fit. Image sizes are 64 X 64 with a pixel size of 1 X 1 arcseconds. Visibility forward fit images shown here were done 
using one elliptical (strong southern footpoint) and circular (weak northern footpoint) Gaussian. 



smaller source sizes, which could be the tendency of th e 
algorith m to over-resolve sour ces (Schm ahl et al J 120071 ). 
PIXON (jPina fc Puetterlll993HMetcalf et a l. 1995) gave 
sour ce sizes similar to t hose of X-ray visibility forward 
fit. IDennis fc Pernakl (J2009D hav e also analysed thi s 
event and confirmed the finding of IKontar et al.l (J2008D . 
We choose to forward fit a circular Gaussian source for 
the northern footpoint and an elliptical Gaussian source 
(Equation [2]) for the southern footpoint to the visibili- 
ties (Equation [TJ, the image shown is a reconstruction of 
the fit results. These fits are shown in Figure [3] and will 
be discussed in detail in i)2.2l Assuming two elliptical 
sources, the weaker source forward fit parameters have 
rather large error bars suggesting that Northern foot- 
point is not sufficiently well-constrained by the data to 
be fitted as an elliptical source. In addition, at the en- 
ergies above 40 keV the weak source is indistinguishable 
from circular. 

The comparison of the images and visibility fit results 
in Figure [5] shows that visibility forward fit gives images 
similar to the ones inferred in other algorithms, although 
there arc differences pointed out above. Despite the dif- 
ferences between the algorithms, all image reconstruc- 
tion algorithms show that the southern footpoint has a) 
a clear elliptical shape, b) the shape of the source be- 
comes more elliptical with growing energy c) the size of 
the source decreases with energy. The northern footpoint 
is also gett ing smaller with ene r gy similar to the souther n 
footpoint (jKontar et al.l 120081 : IDennis fc Pernakl I2009T ). 
but due to lower count rate in the source we cannot reli- 
ably measure the shape of this source. In addition, since 
the northern source is not seen above ~ 100 keV, it could 
be par tially occulted or ha ve stronger magnetic conver- 
gence (Schmahl et al. 2006) with the energetic electrons 
precipitating less to dense layers of the chromosphere, 
producing a fainter footpoint. 

2.2. Characteristic sizes and foot-point locations 

the CLEAN beam size. Note that this correction is only applicable 
for this particular event and cannot be used universally. 



We focus on the brighter southern footpoint fitted with 
an elliptical gaussian (Equation [2]) as more spatial in- 
formation can be accurately recovered compared to the 
northern footpoint. This is a more realistic interpreta- 
tion of the footpoint shape (se e Figure [H compared to 
the previously used circular fit (jAschwanden et aL1l2002t 
IKontar et al.ll2008f) . 

Each forward fit to X-ray visibilities using an elliptical 
gaussian produces 6 parameters given by Equation ([2]): 
positions xo(e) and yo(e) of the X-ray flux maximum (of- 
ten called centroid position), full width half maximum 
FWHM(e), eccentricity e(e) and position angle 6(e), 
along with error values for each of the parameters. Xo(e) 
and yo(e) are measured from disk centre and the posi- 
tion angle 6(e) is the angle between the North-South line 
and the semi-major axis of the ellipse. Multiple sources 
are fitted simultaneously and we fit the weaker north- 
ern source with a circular Gaussian. Both sources can 
be fitted with elliptical sources producing similar results 
for the southern footpoint but highly inaccurate results 
for the northern footpoint. The visibility amplitudes and 
fits as a function of RMC and spacecraft roll angle are 
shown in Figure [3J We used between 6 (course grids) and 
12 (fine grids) visibilities (spatial Fourier components) 
(Figure [3j). The single circular Gaussian fit shows the 
largest amplitudes of normalised residuals. Two circu- 
lar Gaussian fit has smaller amplitudes, but larger than 
the fit using an elliptical and circular Gaussian fits. The 
circular plus elliptical fits adequately reproduce the mea- 
sured photon flux for various roll bins and collimators. 
We note that at the lowest energies 18 — 22 keV, both two 
circular and elliptical plus circular Gaussians give almost 
identical results. Indeed, both footpoints at 18 — 22 keV 
look symmetrical (Figure [2]) . The largest deviations of 
the fits from the data is found in the coarsest grid and 
at the lowest energy (Figure [3]), which could be caused 
by the large scale source (> 36"), probably softer X-ray 
emission from the loop. 

Using the (xo(e),yo(e)) "centroid positions" of the 
source for each energy range, the radial height of the 
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Fig. 3. — Observed X-ray visibility amplitudes (crosses with er- 
ror bars) as a function of subcollimator (2-7) and position angle 
between and 180° of the grids in the stated energy range. The 
red line shows the fitted model using single circular source, the 
green line shows two circular Gaussian fits, and the blue dashed 
line shows an elliptical (southern footpoint) and circular (northern 
footpoint) Gaussian. The bottom panels show normalised residuals 
for the corresponding fits. 



source from the solar centre can be readily determined 
for every energy using: 



R(e) =^o(e) 2 +2/o(e) 2 



(3) 



The semi-major and semi-minor axes of our elliptical 
gaussian fit, a(e), 6(e) respectively are related to the 
FWHM: 

FWHM(e) = v /a(e)6(e), (4) 

and to the eccentricity by: 



(5) 
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Fig. 4. — X-ray visibilities forward fit parameters as a function 
of energy: radial distance of the sources, FWHM of the sources, 
eccentricity, semi-minor and semi-major axes for the northern foot- 
point, and position angle (the angle between semi-major axis and 
North-South direction). 



2.3. Height of X-ray sources above the photosphere 

Let us consider the evolution of the electron flux spec- 
trum in the chromosphere F(E, s) along magnetic field 
lines s using purely collisional transport and ignoring col- 
lective eff ects and effects co nnected with the magnetic 
mirroring ([Brown et al.ll2002l ) . In this a pproximation the 
electron flux spectrum can be written (|BrownllI97lD 



e(e) = v/1 - 6(e)Vo(e) 2 . 
Equations Q and ([5]) can be solved to find 

o(e) = F WHM(e)(l-e(e))~ 1/4 , 6(e) = FWHM(e)(l-e(e 

(6) 
As the analysed flare is right on the Eastern limb and 
close to the solar equator, a(e) and 6(e) correspond to the 
source sizes parallel (width) and perpendicular (vertical 
extent) to the solar surface. The results of the X-ray vis- 
ibility forward fit parameters are summarised in Figure 
01 They again show the trend seen in the images (Fig- 
ure [2]) of decreasing height and source size with energy. 
We can also see that the source becomes more elliptical 
at higher energies, starting with e ~ 0.5 for 18-22 keV 
but increasing to e ~ 0.85 for 75-250 keV. The circu- 
lar Gaussian source fitted to the northern footpoint also 
shows the general trend of decreasing source height and 
FWHM at higher energies but with considerably larger 
errors due to this source being weaker. 



F(E,s)=F (E ] 



E 



(7) 



ere Fq(E) is the injected spectrum of energetic elec- 
trons, taken to be a powerlaw of Fq(E) ex. E~ 5 and 



E (E, 



E 2 + 2K / n(s')ds' 



(8) 



where K = 2we 2 In A, In A is the Coulomb logarithm, e is 
the electron charge. The chromosphere below the tran- 
sition region c an be conveniently assumed to be neutral 



(jBrownllf973l: IKontar et all [20021: ISu et all [2001 the re- 
fore In A = lnA eff = 7 (e.g. IBrowrJll973t lEmslielll978T ). 

The X-ray flux spectrum emitted by the energetic elec- 
trons in a magnetic flux tube of cross-sectional area A 
and observed at 1AU is given as 



I(e,s) 
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F(E,s)a{E,e)dE, (9) 
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Fig. 5. — Upper panel: X-ray spectrum of January 6, 2004 flare 
(solid black histogram) for the accumulation time 06:22:20-06:23:00 
UT (see Figure [TJ and forward fitted spectrum (yellow line). The 
photon spectrum was fitted using a thermal (red dotted line) plus 
thick-target model (blue dashed line). The gray histogram shows 
the background level. Horizontal black and green lines show the 
photon flux spectrum from the southern and northern footpoints as 
found in X-ray visibility forward fits (Figure [2J respectively. The 
lower panel indicate the results of the forward fit of the source 
height with collisional model. The fit results are i?o = 975.2", and 
ho = 155 km. 



where a(E, e) is the isotropic bremsstrahlung cross- 
section, R is the Sun-Earth distance, A is the cross- 
sectional area of the loop. The X-ray flux spectrum 
expressed by Equation @ has a maximum or equiva- 
lently dl(e, s)/ds = for every energy e because of the 
growing density along electron path and si multaneously 
decreasing electron flux due to collisions dBrown et al.l 
l200l lAschwanden et al1l200l . 

Assuming a hydrostatic density profile of 



n(h = r — rn) = no exp 



h 



(10) 



where r is the radial distance from the Sun centre, the 
photospheric density n o = 1.16 x 10 17 cm -3 [fixed value 
(jVernazza et al.l Il981h ] . and rn is the reference height, 
we can find these two free parameters /in and r — rn by 
forward fitting the measured radial distance of maxima 
(Figure [SJ bottom panel) to the model predicted maxima 
by the derivative of equation[9] The height of the sources 
can be found by subtracting the reference height, ro, 
from the radial measurements: 



h(e) = r(e) - r . 



(11) 



To calculate the reference height we assu med the density 
at th e photospheric level to be known (|Vernazza et al.l 
119811) . This helps to remove substantial uncertainties 
related to the reference height of the previous studies (c . f . 
lAschwanden et al.ll200l ILiu et al.ll200l IMrozckj [2006). 

We find a spectral index of S = 3.2 from the spatially 
integrated spectrum, shown in the top panel of Figure [5j 
Forward fitting using a hydrostatic density profile gives 
density scale height of ho = 155 ± 30 km and reference 




10 



Energy [keV] 



100 



Fig. 6. — Top panel: FWHM vertical size of the footpoint source 
(black crosses) as a function of energy; FWHM size given by the 
thick-target model (Equation \$$ (blue dashed lines), multi-thread 
chromosphere (green lines). Bottom panel: height of the maximum 
X-ray emission (black crosses) versus energy. 



height ro = 975.2 ± 0.2". From only fitting a circular 
gaussian to the southern footpoint, instead of elliptical 
to southern with circular to the northern as done in this 
paper, it was previ ously found that h = 140 ± 30 km and 
r = 975.3 ± 0.2" (|Kontar et al.ll2008t) . 

2.4. Vertical extent of the footpoint 

The characteristic size of the source in the vertical di- 
rection (semi-minor axis in our fit) can be straightfor- 
wardly estima ted using the collisional think-target model 
(|Brownlll97l"l) to find the FWHM size from Equation M 
using the density profile found in the previous section. 
Comparing the measured vertical FWHM size and the 
prediction of the length from the thick-target model, we 
see that the measured extent is around 3 — 6" and is 
3 — 6 times larger than the theoretical width (Figure |5J) ■ 
The discrepancy is substantial and cannot be explained 
in terms of the assumed X-rays source model or error 
bars. 

There are a few plausible explanations which can be 
given for the observed vertical FWHM size. One is 
that the structure of the chromosphere may not be uni- 
form over the footpoint cross-section and the footpoint 
is the ensemble of many thin threads (F igure Et as often - 
seen in high reso l ution o ptical images (|Lin et al. 2005; 
iDe Pontieu et al.l l2007t iBerkebile-Stoiser et al.l [2009D . 
Moreover, it is likely that the deposition of electron en- 
ergy can lead to heating of the chromosphere upwards 
with the resulting expans ion changing the density struc- 
ture (e.g. ILiu et aLll2009l ). 

Let us consider a simple chromosphcric model in which 
the hydrostatic density scale height, hx 7 is determined by 
the temperature of the chromosphere hx = kbT / ([im p g) 
where T is the temperature, k\, is the Boltzman con- 
stant, m p is the pro ton mass, /x ^ 1.27 is the mean 
molecular weight (e.g. lAschwanden et al.|[2002l ). and g ~ 
2.74 x 10 4 cm s~ 2 is the solar gravitational acceleration. 
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Fig. 7. — Cross-sectionally uniform (left) and multi-threaded 
(right) chromospheres. The multi-thread chromosphere leads to 
the larger vertical sizes of X-ray sources as observed by RHESSI. 

Thus, the measured density scale height of ~ 155 km im- 
plies an average chromospheric temperature of ~ 6500 K. 
However, the solar atmosphere is not a uniform media 
but inst ead is manifeste d in thin threads (< 0.3)" of fil- 
aments (ILin et alJl2005D . sub-arcsecond dynamic fibrils 
( De Pontieu et al.p 007) and fine structure of microflarcs 
( Berkebile-Stoiser et al.l I2009T ) in the chromosphere with 



different density profiles in each thread (Figure [7]) . Fol- 
lowing the multi-thread model for the magnetic loop in 
the chromosphere, we assume that the energetic electrons 
propagate along different threads with varying density 
profiles and temperatures in the range from 2200 K up 
to 20000 K, corresponding to density scale heights be- 
tween ~ 50 km and ~ 500 km. The average hard X-ray 
flux from many thin threads will be the measured X-ray 
distribution. Averaging X-ray emission from a hundred 
thin threads with temperatures drawn randomly from 
the above range, we successfully reproduce the observed 
vertical sizes (Figure EJ. 

3. DISCUSSION AND CONCLUSIONS 

Using X-ray visibility forward fits, we inferred not only 
the characteristic sizes and positions but the shapes of 
HXR sources. The January 6th 2004 event indicates 
an overall decrease in size of the source and increase of 
the source ellipticity with energy. The FWHM of the 
southern source decreases from ~ 6.5" down to around 
~ 4" while the ellipticity of the source grows from 0.5 
up to 0.9. The source is elongated along the limb as ev- 
ident in nearly zero angle between semi major-axis and 
the limb, such orientation of the source is observed for 
all energy ranges. The vertical extent of the source is 
decreasing by a larger fraction (from 6" down to 3") 
than the horizontal size (from 7" down to 5.5") lead- 
ing to larger elongation of the source along the limb. 
Hence the FWHM of the magnetic flux tube contain- 
ing energetic electrons (semi-major axis) changes from 
~ 9.7 Mm at height ~ 1.1 Mm down to ~ 7.6 Mm at 
0.6 Mm above the photosphere. The northern footpoint 
is fainter, but shows a similar trend: the higher ener- 
gies appear at low heights and the size of the source 
is decreasing with energy suggesting convergence of the 
magnetic field lines along which electrons propagate. Us- 
ing X-ray visibilities we also re-analyzed a flare that 
occurred on February, 20th, 2002 that was previously 
studied by forward fittin g time-modulated lightcurves 
(jAschwanden et all 120021 ). Although this flare is rather 
weak, we found similar results: the higher energy sources 
appear at lower heights. The uncertainties are larger 
than in the Januar y 6th, 2004 event but t he density 
model proposed by lAschwanden et all (|2002D is within 



our error bars. We also note that the size of the sources 
in the February 20th flare decreases with energy similar 
to the event on January 6th. 

Analysing HXR emission from footpoints we found 
that various imaging algorithms (PIXON, Visibility For- 
ward Fit, CLEAN, MEM-NJIT) give generally similar 
results for spatial distributions of X-ray footpoints. Al- 
though RHESSI imaging algorithms can be adjusted by 
the parameter choice and hence X-ray images could be 
somewhat altered, there is a general trend for the algo- 
rithms. CLEAN with default set of parameters has a ten- 
dency to provide larger sources while MEM-NJIT tends 
to over-resolve X-ray sources. PIXON and visibility for- 
ward fit show very similar results. Visibility forward fit 
allows us to study sub-arcsecond distribution of hard X- 
ray sources in suitably orientated bright flares. Due to 
systematic differences in sizes we obtained for January 
2004 and February, 2002 flares, we suggest that CLEAN 
and MEM-NJIT should be used for source size/shape 
measurements of X-ray sources with extreme caution. 

The northern footpoint could be partially occulted as 
the highest energy photons come predominantly from the 
southern footpoint, but we note that this bright foot- 
point is unlikely to be occulted. As pointed out by G. 
HurforcQ, a source partially subtended by the solar disk 
should have a sharp edge, where the brightness of the 
source will drop from maximum to zero over rather small 
radial distance. The derivative of the source brightness in 
the x-dircction (perpendicular to the limb), dl{x,y)/dx 
will have a maximum at the limb, where I(x, y) has a 
sharp drop. It is evident from Equation ([1]) that the cor- 
responding visibilities, 2niuV(u,v), should have a well 
pronounced maximum, which should be evident in the 
measured amplitudes of visibilities. Specifically, the finer 
grid RMCs should show large amplitudes when the grids 
are parallel to limb, i.e. the visibility amplitudes are 
much larger at the phase angles 0° and 180°, which 
is not evident in the event under study (cf Figure [3]). 
To make the discussion more complete, we note that the 
line of sight effects cannot be definitively ruled out. The 
footpoints might be projections of two rather long flare 
ribbons viewed almost parallel to the line of sight, so that 
the vertical extension is the projection of different height. 
Finally we should note that if the occultation height is 
small, < 0.5", our height measurements are lower lim- 
its, but the major conclusion about the vertical extend 
is the same. If, though unlikely, the lower part of the 
loop (footpoints) is occulted, these observations provide 
an interesting question to the flare models as to why the 
sources sizes decrease with energy and the higher energy 
sources appear lower and not at the lowest visible loca- 
tion where the density is the highest. 

Our measurements show that while the locations of 
the maxima of X-ray emission are consistent with simple 
collisional transport in single density scale height chro- 
mosphere, the vertical sizes do not agree with the as- 
sumption of field aligned electron transport. The ver- 
tical extent of X-ray sources is 3-6 times larger than in 
the purely collisional model in single-density-scale-height 
chromosphere. However, a chromospheric model involv- 
ing multiple density threads within the flux tube of a 

2 P resentati on at 9th RHESSI workshop in Genoa 

|http://sprg. ssl.berkeley.edu/~krucker/gcnoa/position/XrayLimb-Genoa.ppt | 



The sub-arcsccond hard X-ray structure of loop footpoints in a solar flare 



footpoint can explain both the position of the maximum 
and the vertical size of the sources. We note that pitch 
angle scattering due to Coulomb collisions is likely to be 
insufficient to produce so strong expansion. The X-ray 
source size increase due to collisional pitch angle scatter- 
ing will be about a quarter of the e lectron stopin g depth 
for initially field aligned electrons (|Conwavll200ClD . How- 
ever, strong non- collisional scattering or wave-particle 
interactions (e.g. lHannah et al.l [2009? ) might boost the 
vertical source sizes to the measured value and hence 
cannot be excluded and will also be consistent with lack 
of downward anisotropy found in X-ray flare emission 
(|Kontar fc Brownll2006l : lKasparova et aLJl2007D . We note 
that the adopted model does not account for the mag- 



netic field and its effects on particle transport, which 
could lead to larger source sizes and is subject of addi- 
tional modeling. 
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